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Quantum  transport:  Silicon  inversion  layers  and  InAIAs-InGaAs 
heterostructures 

D.  Vasileska,  T.  Eldridge,^'  and  D.  K.  Ferry 

Center  for  Solid  State  Electronics  Research,  Arizona  State  University,  Tempe,  Arizona  85287-(>2U0 
(Received  21  January  1996;  accepted  26  April  1996) 

Using  a  real-time  Green’s  functions  formalism,  we  investigate  the  influence  of  depletion  charge 
scattering  on  the  room  temperature  mobility  in  scaled  silicon  metal-oxide-semiconductor  field-effect 
transistors  and  low-temperature  transport  in  Ino.4Alo.6As-Ino4Gao.6As  modulation  doped 
heterostructures.  Our  simulation  results  suggest  that  depletion  charge  scattering,  which  is  usually 
ignored,  has  considerable  impact  on  the  electron  transport  in  silicon  inversion  layers  near  the 
threshold  gate  voltage,  even  at  room  temperature.  We  also  find  that  the  weighting  coefficients  a  and 
b  (for  the  inversion  and  depletion  charge  densities)  strongly  depend  on  the  substrate  doping  and 
deviate  from  that  reported  in  the  literature.  In  the  case  of  modulation  doped  heterostructures,  the 
low-temperature  mobility  is  limited  by  alloy  and  Coulomb  scattering.  Intersubband  scattering 
considerably  affects  the  broadening  of  the  states  which,  in  turn,  leads  to  mobility  reduction. 
©  1996  American  Vacuum  Society. 


I.  INTRODUCTION 

Electron  transport  properties  of  quasi-two-dimensional 
(Q2D)  systems  have  been  studied  quite  extensively,  experi¬ 
mentally  as  well  as  theoretically.  The  inversion  layer  at  a 
silicon  interface'  represents  the  most  widely  employed 
physical  realization  of  the  Q2D  concepts  and  the  electron 
mobility  of  the  inversion  layer  electrons  is  one  of  the  most 
important  parameters  for  device  characterization.  It  has  been 
confirmed  experimentally^'^  that  for  substrates  with  low  dop¬ 
ing  densities,  the  electron  mobility  follows  a  universal  curve 
when  plotted  as  a  function  of  the  effective  electric  field. 
However,  pushing  the  devices  towards  scaling  limits''  neces¬ 
sarily  leads  to  thinner  oxide  films  (<10  nm)  and  higher  sub¬ 
strate  dopings  (>10'’  cm"^).  For  these  devices,  it  was  re¬ 
cently  demonstrated^  that  experimental  mobilities 
significantly  deviate  from  the  universal  curves  near  the 
threshold  voltage  as  a  consequence  of  significant  depletion 
charge  scattering. 

Since  the  invention  of  the  modulation  doping  technique, 
transport  properties  of  two-dimensional  electron  gases  in 
high-mobility  modulation  doped  AlGaAs-GaAs  heterojunc¬ 
tions  (MDHs)  have  also  been  a  subject  of  great  interest  from 
both  the  technological  and  fundamental  point  of  view.  “ 
The  enhanced  mobilities  in  these  structures  are  attributed  to 
the  spatial  separation  of  the  electrons  from  their  parent  do¬ 
nors,  which  significantly  reduces  the  scattering  rate  from  the 
remote  impurities."”'^  This  has  stimulated  much  research 
toward  their  application  in  low-power,  high-speed  devices. 
Recently,  there  has  been  an  increased  interest  devoted  to  the 
growth  of  InAlAs-InGaAs  lattice  matched 
heterojunctions'''"'®  which  have  better  carrier  confinement 
due  to  a  larger  conduction-band  discontinuity  and  a  lighter 
eonduction-band  effective  mass  (compared  to  AlGaAs- 
GaAs).  However,  since  the  electron  gas  in  the  InAlAs- 
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InGaAs  MDH  is  mostly  localized  within  a  ternary  random 
alloy  (InGaAs),  alloy  scattering^'^  is  much  more  effective  so 
that  it  is  hard  to  imagine  that  the  low-temperature  mobility  of 
InGaAs  based  modulation  doped  heterostructures  may  ex¬ 
ceed  the  very  high  mobilities  already  reported  in  AlGaAs- 
GaAs  heterostructures. 

The  theoretical  model  for  the  calculation  of  the  subband 
structure  and  the  electron  mobilities  of  the  Q2D  electron  gas, 
taking  into  consideration  all  major  scattering  mechanisms 
and  both  intrasubband  and  intersubband  scattering,  is  briefly 
described  in  Sec.  11.  In  Sec.  Ill,  devoted  to  transport  in  sili¬ 
con  inversion  layers,  we  describe  the  influence  of  depletion 
charge  scattering  on  the  mobility  results  near  the  threshold 
gate  voltage.  We  also  show  how  the  weighting  coefficients 
for  the  inversion  and  depletion  charge  densities,  which  ap¬ 
pear  in  the  definition  of  the  average  electric  field,^^  vary  with 
doping  density.  Low-temperature  transport  properties  of  the 
modulation  doped  Ino.4Alo.6As-Ino.4Gao.6As  heterostructures 
are  discussed  in  Sec.  IV.  Some  concluding  comments  related 
to  the  present  work  are  given  in  Sec,  V. 


II.  THEORETICAL  MODEL 

The  subband  structure  is  obtained  from  a  self-consistent 
solution  of  the  Schrodinger  and  Poisson  equations,  coupled 
with  the  Dyson  equation  for  the  retarded  Green’s  function. 
Exchange-correlation  corrections  to  the  chemical  potential 
are  calculated  within  the  local  density  approximation.  The 
self-energy  contributions  from  the  various  scattering  mecha¬ 
nisms  are  evaluated  by  employing  the  self-consistent  Bom 
approximation.  Screening  is  treated  within  the  random  phase 
approximation.  In  the  calculation  of  the  conductivity,  we 
have  taken  account  of  a  correction  term  due  to  the  particle- 
hole  ladder  diagrams.  More  details  on  the  calculational 
model  and  the  analytical  expressions  for  the  broadening  of 
the  states  and  conductivity  can  be  found  elsewhere. 
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Fig.  1.  Room  temperature  mobility  vs  inversion  charge  density  for  two 
different  devices  with  high  doping  densities.  For  both  devices,  interface-trap 
density  is  A^it=1.75Xl0‘^  cm~^.  We  use  an  exponential  model  for  the 
surface-roughness  autocorrelation  function.  We  find  that  the  rms  height  of 
the  bumps  is  0.26  nm  (0.25  nm)  for  the  device  with  lower  (higher)  doping. 
For  both  devices,  roughness  correlation  length  is  estimated  to  be  1.4  nm. 


III.  SILICON  INVERSION  LAYERS  WITH  HIGH 
DOPING  PROFILES 

The  physical  model  employed  in  these  calculations  in¬ 
cludes  scattering  with  all  bulk  phonons  (acoustic  and  nonpo¬ 
lar  optical),  surface-roughness  scattering,  and  Coulomb  scat¬ 
tering  from  both  ionized  impurities  in  the  depletion  region 
and  interface/oxide  charges.  For  intervalley  scattering  (zero- 
order  and  first-order  interactions),  we  use  the  same  phonon 
energies  and  coupling  constants  as  those  relevant  for  bulk 
silicon.^®  The  anisotropy  of  the  deformation  potential  inter¬ 
action  is  also  accounted  for.^^  More  details  about  the  model 
can  be  found  in  the  Appendix. 

In  Fig.  1,  we  compare  our  simulation  results  for  the  mo¬ 
bility  with  the  experimental  effective  mobility  data  for 
rt -channel  metal-oxide-semiconductor  field-effect  transistors 
(MOSFETs)  fabricated  on  (100)  Si  wafers.  Filled  circles  and 
triangles  represent  the  experimental  effective  mobility  data 
derived  from  drain  conductance  measurements.^  The  inver¬ 
sion  charge  density  has  been  determined  from  the  gate- 
channel  capacitance  measurements.  Solid  (dashed)  lines  are 
the  corresponding  simulation  results  when  depletion  charge 
scattering  is  included  (omitted)  in  the  calculations.  A  total  of 
four  subbands  (three  from  the  unprimed  and  one  from  the 
primed  ladder  of  subbands)  are  taken  into  account.  We  use 
the  real  density  of  states  function  (DOS)  when  evaluating  the 
self-consistent  potential-energy  profile.  In  contrast  to  some 
previous  studies,  we  do  not  ignore  the  real  self-energy 
corrections  (renormalizations  of  the  spectrum)  when  evaluat¬ 
ing  the  broadening  of  the  states  and  conductivity.  We  find 
that  depletion  charge  scattering  leads  to  mobility  reduction  of 
about  9%  near  the  gate  voltage  thresholds  (small  values  of 
N^)  for  the  device  with  A^^=3X10^^  cm"^.  For  the  device 
with  higher  substrate  doping  (A^^  =7.7X10^^  cm”^),  the  mo¬ 
bility  reduction  is  about  12.5%.  For  both  devices,  the  impor¬ 
tance  of  surface-roughness  scattering  is  shifted  towards 
lower  gate  voltages.  In  addition,  scattering  due  to  acoustic 


Substrate  doping  [cm"^ 


Fig.  2.  Variation  of  the  weighting  coefficients  a  and  b=\-{z)/W  with 
substrate  doping.  We  find  that  the  penetration  of  the  wave  functions  into  the 
oxide  is  about  1  nm  for  very  high  values  of  .  Slightly  higher  values  for 
the  coefficient  b  were  obtained  for  V,>1X10^^  cm"^  (not  shown  in  this 
figure)  due  to  the  closer  confinement  of  the  carriers  to  the  interface  and 
increased  population  of  the  lower-lying  subbands. 

and  nonpolar  optical  phonons  is  found  to  be  of  lesser  impor¬ 
tance  in  these  structures. 

It  has  already  been  reported^^’^^  that  room  temperature 
electron  mobilities  show  universal  behavior  independent  of 
the  substrate  doping  when  plotted  as  a  function  of  the  aver¬ 
age  electric  field  E^^={e/e^^)(aNs-\-bN^^pi),  where  a  and  b 
are  the  weighting  coefficients  for  the  inversion  and  deple¬ 
tion  A^depi  charge  densities.  In  Fig.  2  we  show  the  variation  of 
these  two  coefficients  with  substrate  doping.  We  find  that  a 
significantly  deviates  from  0.5  for  substrates  with  very  high 
doping  densities.  We  also  observe  that  the  weighting  coeffi¬ 
cient  b  decreases  with  increasing  substrate  doping  due  to  a 
decrease  of  the  depletion  region  width.  Results  shown  in  Fig. 

2  differ  from  those  obtained  by  Matsumoto  and  Uemura^"^ 
(a  =0.5  and  ^  =  1),  for  example,  because  in  their  analytical 
result  for  the  effective  electric  field,  the  terms  of  the  order  of 
{z)IW  (suggested  by  Krutsick  and  White^^),  where  (z)  is  the 
centroid  of  the  electron  density  distribution  and  W  is  the 
depletion  region  width,  were  neglected. 

IV.  lno.4Alo6As-lno4Gao6As  MODULATION  DOPED 
HETEROSTRUCTURES 

At  low  temperatures,  the  dominant  scattering  mechanisms 
in  the  Ino.4Alo6As-Ino.4Gao.6As  MDH  shown  in  Fig.  3  are: 
remote  impurity  (from  the  ^depletion  layer),  interface-trap, 
alloy  disorder,  interface-roughness  and  deformation  poten¬ 
tial,  and  piezoelectric-coupled  acoustic  mode  scattering. 
Polar  optical  phonon  scattering  is  found  not  to  play  any  sig¬ 
nificant  role  at  1.6  K.  Screening  of  the  alloy  disorder^^  and 
deformation  potential  scattering  potentials  by  the  free  earn¬ 
ers  is  also  neglected  because  of  their  short-range  nature.  This 
effect  should  certainly  have  been  included  for  the  long-range 
electron-piezoelectric  acoustic  phonon  interaction,  but  since 
this  scattering  process  is  extremely  weak  at  1.6  K,  its  matrix 
element  is  also  left  unscreened.  Assuming  that  deformation 
potential  scattering  is  an  elastic  process,  one  has  ho)Q/k^T 
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Gate 

30  nm  ud-Ino.4Alo.6As  (barrier  layer) 
10  nm  ud-Ino.4Alo.6As  (spacer) 

50  nm  ud-Ino.4Gao.6As  (well) 

0.1 /an  ud-Ino.4Alo.6As  (buffer  layer) 
0.1 /rm  ud-Ino.3Alo.7As  (buffer  layer) 
0.1 /im  ud-Ino.2Alo.8As  (buffer  layer) 
0.1 /im  ud-Ino.1Alo.9As  (buffer  layer) 
0.1  /rm  ud-  GaAs  substrate 
Semi- insulating  GaAs  substrate 


Fig.  3.  Cross-sectional  view  of  the  <5-doped  Ino.4Alo.6As-Ino.4Gao.6As  hetero- 
structure  on  (lOO)-oriented  semi-insulating  GaAs  substrate  with  intermedi¬ 
ate,  undoped,  and  compos itionally  step-graded  buffer  layers. 


--  yl3c^m^/pksf  =  where  o)q  is  the  phonon  fre¬ 

quency,  kg  is  the  Boltzmann  constant,  is  the  longitudinal 
elastic  constant,  m*  is  the  effective  mass,  and  p  is  the  crystal 
density.  Using  the  material  parameters  summarized  in  Table 
I,  we  find  that  the  equivalent  temperature  is  T^^O.22  K. 
Thus,  the  equipartition  approximation  can  be  assumed  to 
hold  for  temperatures  above  1  K  for  the  acoustic  mode  scat¬ 
tering  of  thermal  electrons.  The  matrix  elements  for  surface- 
roughness,  acoustic  phonon,  and  alloy  disorder  scattering  can 
be  found  elsewhere.^^  For  the  wavevector  dependent  matrix 
element  for  scattering  between  subband  n  and  m  due  to  Cou¬ 
lomb  charges  located  at  a  distance  Zi<0  from  the  interface 
(in  the  InAlAs  layer),  we  use  the  following  expression: 


Table  I.  Values  of  the  In^Gai  _;,As  material  parameters  used  in  the  calcula¬ 
tion  of  the  broadening  of  the  states,  density  of  states  function,  and  mobility. 
For  the  electron  effective  mass  and  static  dielectric  constants  of  the 
Ino.4Alo.6As,  we  use  m*=0.99  mo  and  e=11.5  60,  respectively.  The  energy 
gaps  of  strained  and  unstrained  Inj^Gai  -j;As  layers  can  be  found  in  Ref.  25. 


Parameter 

Notation  and 
units 

Value 

Relative  electron  effective  mass 

m* 

0.058  for  a: =0.4 

Relative  static  dielectric  constant 

6 

13.1  +  1.5a: 

Deformation  potential  constant 

D  (eV) 

8.6-2.8X 

Piezoelectric  constant 

P 

0.052-0.007  9x 

Longitudinal  elastic  constant 

(XIO®  N/cm^) 

13.97-3.99x 

Lattice  constant 

a  (A) 

5.653  25+0.405  15x 

Alloy  disorder  parameter 

(V>  (eV) 

0.5-0.55 

Crystal  density 

p  (gW) 

5.504 
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Gate  voltage  [V] 


Fig.  4.  Gate-voltage  dependence  of  the  total  electron  density  for  the  MDH 
shown  in  Fig.  3.  Also  shown  in  this  figure  are  the  experimental  data  derived 
from  the  low-field  Hall  and  resistivity  measurements,  Shubnikov-de  Haas 
magnetoresistance  measurements  as  well  as  the  theoretical  simulations  re¬ 
ported  in  Ref.  16.  The  effective  mass  enhancement  due  to  conduction  band 
nonparabolicity  for  the  carriers  residing  in  the  ground  and  first  excited  sub¬ 
band  is  shown  in  the  inset. 


{n\U{q,Zi)\m)=^^  j _Jz  Uz) 


+  q\z+ii\  (^^(2)  + 


dz  <pn{z)e 

J  -00 


q(€i,+  €j 


in  which  the  image  term  due  to  the  difference  of  the  static 
dielectric  constants  in  the  well  and  in  the  barrier  layer 
(£(,)  is  properly  taken  into  account.  The  effective  mass  en¬ 
hancement  with  increasing  gate  voltage  due  to  conduction 
band  nonparabolicity  is  also  taken  into  account. 

When  solving  for  the  subband  structure,  the  number  of  the 
self-consistent  field  iterations  was  considerably  decreased 
when  using  the  analytical  solutions  of  the  Poisson  equation. 
We  calculate  that  in  equilibrium,  the  electron  transfer  from 
the  <5-doped  region  to  the  well  (see  Fig.  3)  is  governed  by 

EF+^o-eVc  =  ^E,+  — - —  zif^i(z)dz 

-  — (W+Z,)2  l^r  1-  f  fl'^(z)dz 

€b  l  i  [  J-oo 


(2) 


where  L=30  nm  and  1V  =  10  nm  are  the  widths  of  the  barrier 
and  spacer  layers,  Ep  is  the  Fermi  level,  4>o=0.8  eV  is  the 
Shottky  barrier  height,  Vc  is  the  gate  voltage,  AEj.=0.52  eV 
is  the  conduction  band  offset,  A^£,=3X10*^  cm  ^  is  the  Si 
5-doping  sheet,  and  is  the  sheet  electron  density  of 

the  ith  subband  in  the  barrier  (quantum  well). 
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Fig.  5.  Cumulative  and  component  theoretical  mobilities  vs  gate  voltage. 
Filled  circles  are  the  experimental  data  determined  from  Hall  effect  and 
resistivity  measurements  (Ref.  16).  The  material  parameters  employed  in  the 
simulations  are  summarized  in  Table  I.  Also  shown  in  the  figure  are  the 
mobilities  of  the  carriers  in  the  ground  (ytti)  and  first  excited  subband 


Self-consistent  results  for  the  gate-voltage  dependent 
electron  densities  are  shown  in  Fig.  4.  We  observe  good 
agreement  between  the  experimental  data  derived  from  the 
Hall  effect  and  SdH  measurements^^  and  the  theoretical  cal¬ 
culations  for  all  gate  voltages.  Our  simulation  results  are  also 
in  agreement  with  the  theoretical  calculations  reported  by 
Wieder  et  al.  We  find  that  the  significant  effective  mass 
enhancement  (see  the  inset  of  Fig.  4)  due  to  conduction  band 
nonparabolicity  prevents  the  population  of  the  third  subband 
until  Vq>0A  V.  We  also  find  that,  in  contrast  to  silicon  in¬ 
version  layers  where  exchange-correlation  effects  are  signifi¬ 
cant  even  at  room  temperature,  many-body  corrections  do 
not  considerably  affect  the  subband  structure  in  the  MDH 
because  of  the  small  effective  mass. 

Gate- voltage  dependence  of  the  low-temperature  theoreti¬ 
cal  mobilities  is  shown  in  Fig.  5.  Although  the  average  mo¬ 
bility  is  ifii/N s ,  measured  mobilities  shown  in  Fig.  5 

are  Hall  mobilities'^  so  that  for  a  multisubband  case,  the 
expression  is  more  suitable  for  compar¬ 

ing  the  experimental  results  with  the  theoretical  calculations. 
To  take  into  account  the  diffusion  effects,  the  ^doping  sheet 
is  modeled  as  a  4  nm  wide  layer  centered  8  nm  away  from 
the  interface  with  Si  doping  density  of  7.52X10^*  cm“^. 
Interface-roughness  scattering  is  assumed  to  exist  at  both 
heterointerfaces  of  the  quantum  well.  We  assume  that  the 
rms  height  of  the  roughness  is  one  monolayer  (ML)  (A=2.83 
A).  The  lateral  spread  of  the  interface  roughness  is  taken  to 
be  ^=100  A.  For  this  particular  heterostructure,  we  find  that 
surface-roughness  and  acoustic  phonon  limited  component 
mobilities  (not  shown  on  the  figure)  are  of  the  order  of 
(0.5-1)  X 10^  and  10^  cm^A^  s,  respectively.  By  fitting  the 
calculated  mobilities  at  high  gate  voltages  to  the  experimen¬ 
tal  data,  we  were  able  to  determine  the  alloy-disorder  scat¬ 
tering  parameter  in  the  system  to  be  in  the  range  0.5-0.55. 
Component  mobility  results  presented  in  Fig.  5  suggest  that 
the  inclusion  of  alloy  scattering,  in  addition  to  Coulomb  scat¬ 
tering  (from  the  ionized  donors  in  the  ^doped  region  and 


interface  traps),  leads  to  good  agreement  between  the  theo¬ 
retical  calculations  and  the  experimental  data.  The  dip  in  the 
theoretical  mobility  at  yG  =  “0-85  V  signals  that  a  new  chan¬ 
nel  (second  subband)  for  scattering  opens  in  the  MDH  which 
leads  to  a  sudden  increase  in  phase  space  for  scattering  be¬ 
tween  the  two  subbands  to  occur.  Since  the  carriers  residing 
in  the  first  excited  subband  are  at  about  four  times  larger 
average  distance  from  the  heterointerface  (compared  to  those 
in  the  ground  subband),  they  do  not  feel  as  much  the  Cou¬ 
lomb  scattering  potentials  and  exhibit  about  two  times  higher 
mobility. 

V.  CONCLUSIONS 

Our  results  suggest  that  depletion  charge  scattering,  which 
has  been  ignored  in  some  previous  theoretical  studies,  is  re¬ 
sponsible  for  the  deviation  of  the  mobility  from  the  universal 
curves  near  the  threshold  voltage.  At  high  gate  voltages  (high 
values  of  due  to  the  closer  confinement  of  the  carriers  to 
the  interface,  the  magnitude  of  the  mobility  in  these  highly 
doped  samples  is  considerably  reduced  by  surface-roughness 
scattering,  in  addition  to  the  still  significant  intervalley  scat¬ 
tering.  We  also  find  that  the  weighting  coefficients  a  and 
b — that  appear  in  the  definition  of  the  average  electric 
field — show  a  pronounced  doping  dependence.  For  the 
InQ  4Alo  6As-Ino  4Gao  6As  MDH,  we  find  that  at  r=1.6  K  the 
magnitude  of  the  mobility  is  limited  by  Coulomb  and  alloy 
scattering.  From  the  subband  structure  (gate-voltage  depen¬ 
dence  of  the  electron  density)  and  mobility  simulations,  we 
calculate  that  the  relevant  effective  mass  at  the  conduction 
band  minima  of  the  Ino.4Gao  6As  ternary  alloy  is  w* =0.058 
niQ.  An  observed  effective  mass  enhancement  suggests  that 
nonparabolicity  effects  are  considerable  in  this  MDH  and 
must  be  taken  into  account  if  accurate  mobility  results  are 
desired.  The  effect  of  the  many-body  interactions,  calculated 
in  the  local  density-functional  approximation,  is  found  not  to 
be  important  for  the  MDH, 
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APPENDIX 

The  expression  for  the  matrix  elements  for  scattering  be¬ 
tween  subbands  n  and  m  due  to  a  Coulomb  charge  in  the 
depletion  region,  in  the  oxide,  or  right  at  the  interface,  in  the 
presence  of  a  different  dielectric  medium,  can  be  found  in 
the  review  article  by  Ando,  Fowler,  and  Stem.^  For  the 
surface-roughness  matrix  element,  we  use  the  result  obtained 
by  Ando^^  who  corrected  the  result  obtained  by  Matsumoto 
and  Uemura^^  by  taking  into  account  the  change  in  the  elec¬ 
tron  density  distribution  and  the  effective  dipole  moment  of 


JVST  B  -  Microelectronics  and  Nanometer  Structures 


2784 


Vasileska,  Eldridge,  and  Ferry:  Quantum  transport:  Si  inversion  layers 


2784 


the  deformed  Si-Si02  surface  in  addition  to  the  usual  local 
potential  energy  term  that  equals  the  linear  term  of  the  Taylor 
expansion  of  the  surface  potential. 

For  anisotropic  elastic  continuum  such  as  silicon,  for 
which  the  deformation  potential  constant  a  is  a  tensor, 
matrix  element  squared  for  scattering  between  subbands  n 
and  m  due  to  acoustic  phonons,  in  the  usual  equipartition 
approximation,  equals  to^^ 

|(«|C/r(q)|m)P=^  (Al) 

where 

[  dz  (A2) 

Jo 

is  the  sound  velocity,  and  p  is  the  Si  mass  density.  The 
effective  deformation  potential  constant  that  appears  in  ex¬ 
pression  (Al)  is  calculated  from 

1 

=  —  dq,  (A3) 

^nm  Jo 


where 

■^nm{dz)=  f  dz  tlfniz)e"^^’^ll/„{z)  (A4) 

Jo 

and  Oq  is  the  angle  between  the  wave  vector  Q  of  the  emitted 
(absorbed)  phonon  and  the  longitudinal  axis  of  the  valley. 
According  to  Herring  and  Vogt,^^  the  effective  deformation 
potentials  A;^(^q)  (X=LA  or  TA)  that  appear  in  Eq.  (A3)  are 
given  by 

ALA(^Q)““Hd+H„  cos^C^q)  (A5a) 

and 

A  ta(  ^q)  H  „  cos(  eq)  sin(  (9q)  .  (A5b) 

In  Eqs.  (A5a)  and  (A5b),  H„  is  the  uniaxial  shear  potential, 
and  is  the  dilatation  potential  that  are  believed  to  have 
values  of  approximately  9.0  and  —11.7  eV,  respectively.^ 

It  is  important  to  point  out  that  Eq.  (A5b)  accounts  for  the 
contribution  of  both  TA  branches. 

The  matrix  element  for  nonpolar  optical  phonon  scatter¬ 
ing  is  generally  found  from  a  deformable  ion  model.  The 
zero-order  process,  which  describes  a  perturbation  that  can 
be  associated  with  the  local  dilatation  or  compression  of  the 
lattice  which  produces  local  fluctuation  in  the  electron  en¬ 
ergy,  dominates  at  moderately  high  electron  energies.  The 
matrix  elements  for  scattering  between  subbands  n  and  m 
(that  belong  to  a  and  ^  valleys)  due  to  this  zero-order  inter¬ 
action  is  given  by 

=  (A6) 

where  is  the  deformation  field  and  Wox  is  the  frequency  of 
the  phonon  mode. 

When  the  zero-order  matrix  element  for  the  optical  or 
intervalley  interaction  vanishes,  then  is  identically  zero. 
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In  this  case,  one  has  to  consider  the  first-order  term  of  the 
interaction.  The  squared  matrix  element  for  scattering  be¬ 
tween  subbands  n  (a  valley)  and  m  (^  valley)  due  to  this 
first-order  process  is  given  by^^ 

where 


0 


(A8) 


In  the  scattering  among  the  equivalent  valleys,  there  are 
two  types  of  phonons  that  might  be  involved  in  the  process. 
The  first  type,  the  so-called  g  phonon  couples  the  two  v^leys 
along  opposite  ends  of  the  same  axis,  i.e.,  (100)  to  (100). 
This  is  an  umklapp  process  and  has  a  net  phonon  wavevector 
O.STr/a.  The  so-called  /  phonons  couple  the  (100)  valley  with 
(010),  (001),  etc.  Hence,  the  reciprocal  lattice  vector  in¬ 
volved  in  the  g -process  is  Gioo  and  that  for  an  /  process  is 
^111- 

Within  a  four  subband  approximation,  scattering  between 
the  two  valleys  in  the  €q,  and  62  subbands  involves  only 
g-type  phonons.  The  scattering  between  these  three  minima 
is  treated  by  using  a  high-energy  phonon  of  750  K  activation 
temperature  (treated  as  zero-order  interaction)  and  134  K 
phonon  treated  via  first  order  interaction.  Scattering  between 
the  60,  61 ,  and  62  subbands  and  the  four  60  subbands  involves 
/  phonons  with  activation  temperatures  of  630  and  230  K 
treated  via  zero-order  and  first-order  interaction,  respectively. 
Scattering  between  the  €q  subbands  involves  both  g  and  / 
phonons  with  activation  temperatures  of  630  K  (zero-order 
interaction)  and  190  K  (first-order  interaction).  All  of  the 
high-energy  phonons  are  assumed  to  be  coupled  with  a  value 
of  Dx-9X10^  eV/cm,  and  all  of  the  first-order  coupled 
phonons  are  assumed  to  be  coupled  with  Z)ix=5.6  eV.  These 
values  are  consistent  with  the  results  given  in  Refs.  20  and 
32. 
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